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The Schrodinger wave functional ip — exp — S{A1(x)} for the d — 3+1 QCD vacuum is a partition 
function constructed in d = 4; the exponent 25* [in \ip\ 2 — exp(— 2S)] plays the role of a d — 3 
Euclidean action. We start from a simple conjecture for 5 based on dynamical generation of a gluon 
mass M in d = 4, then use earlier techniques of the author to extend (in principle) the conjectured 
form to full non-Abelian gauge invariance. We argue that the exact leading term, of O(M), in an 
' expansion of S in inverse powers of M is a d — 3 gauge-invariant mass term (gauged non-linear sigma 

model); the next leading term, of 0(1/M), is a conventional Yang-Mills action. The d = 3 action 
, that is (twice) the sum of these two terms has center vortices as classical solutions. The d — 3 gluon 

mass 7713, which we constrain to be the same as M, and d = 3 coupling g 2 are related through the 
conjecture to the d — 4 coupling strength, but at the same time the dimensionless ratio mg/g 2 can 
be estimated from d = 3 dynamics. This allows us to estimate the d — 4 coupling a s (M 2 ) in terms 
of the strictly d = 3 ratio mzj g\\ we find a value of about 0.4, in good agreement with an earlier 
theoretical value but somewhat low compared to the QCD phenomenological value of 0.7 ±0.3. The 
wave functional for d — 2 + 1 QCD has an exponent that is a d = 2 infrared-effective action having 
. both the gauge-invariant mass term and the field strength squared term, and so differs from the 

conventional QCD action in two dimensions, which has no mass term. This conventional d = 2 
QCD would lead in d — 3 to confinement of all color-group representations. But with the mass 
term (again leading to center vortices), only AT-ality ^ mod N representations can be confined 
©; (for gauge group SU{N)), as expected. 

PACS numbers: ll.15.-q, 12.38.-t, 11.15.Tk UCLA/06/TEP/29 

o : 

^ ■ I. INTRODUCTION 

The functional Schrodinger equation (FSE) for gauge theories, while no simpler to solve (and perhaps harder, in 
some ways) than any other non-perturbative formulation of QCD, has often been used over the years to gain insight into 
various aspects of QCD or, more generally, SU(N) gauge theory 0, 0, i, i, i, 0, 0, 0, 0, El M, G3, EMI El El E3 • 
. [ However, few of these works address the important question of how confinement is expressed in the FSE. 

In any approach to the FSE for QCD that purports to reveal confinement, there are two important prerequisites: 
The first is gauge invariance, and it has been addressed many ways. The second is the need to insure that there 
are only short-range field-strength correlations; otherwise (see, e. g., the qualitative and in some ways incomplete 
discussion of Feynman [6|) there cannot be confinement. Given these, confinement further requires long-range pure- 
gauge contributions to the potential. These long-range pure-gauge parts appear in the FSE as massless longitudinally- 
coupled scalars that mimic Goldstone fields, although of course there is no symmetry breaking in QCD. Just as with 
conventional Goldstone fields, these massless poles do not appear in the QCD S-matrix; this would be so even if QCD 
were not a confining theory. As is well-known, center vortices, solitons of an infrared-effective action for QCD that 
encapsulates dynamical and gauge-invariant generation [H, El] of a gluon mass M, show just these properties and 
so provide a confinement mechanism. This mass has been estimated theoretically from phenomenology [20l. l2l|. 
and on the lattice [13], all yielding values of 600±200 MeV. The center vortices in d = 3 are characterized by closed 
strings that (generically) constitute the constant-time cross-sections of d = 4 center vortices; a confining condensate of 
center vortices in d = 4 is therefore mirrored by a similar condensate in d = 3. (Of course, the classical local minimum 
describing a single or a few center vortices is not relevant in isolation; it is necessary that there be a entropy-driven 
condensate of vortices. We do not discuss that issue here.) 

In d = 3 + 1 the FSE describes four-dimensional dynamics in d = 3 terms, because the exponent S in the vacuum 
wave functional ip — exp(— S) is (half of) a d = 3 action, which we label Id=3, depending on the spatial gauge 
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potentials at zero time. Many authors have discussed center vortices for QCD in strictly d = 4 terms. Our question 
is, how are such solitons — and hence confinement — described in the FSE action 2S — 

Our answer proceeds in four steps. The FSE exponent S is an infinite series of n-point functions integrated over 
the spatial components of n gauge potentials (see the Appendix, which reviews earlier work [8[ on the FSE, as well as 
Sec. IIII[) . The first step, described in Sec. [HI considers the lowest-order term S2 of this expansion, which is quadratic 
and shows only Abelian U(1) N _1 gauge invariance. Our conjectured form of S2 exactly satisfies the FSE with an 
Abelian gauge Hamiltonian that phenomenologically describes a gauge-invariant gluon mass M; it is essentially N 2 — 1 
copies of the Abelian Higgs model with infinite Higgs mass. 

Since our focus is on confinement, an infrared phenomenon, we will use techniques and approximations that are useful 
in the infrared regime, even though they may misstate ultraviolet-dominated phenomena. In particular, although we 
treat the gluon mass M as a constant, it is actually a running mass M(k 2 ) evaluated on-shell. In order that there 
be dynamical mass generation in QCD, the running mass must vanish for large momentum k 2 [Hj]. This vanishing 
cures certain short-distance singularities of the center-vortex solitons coming from an infrared-effective action. We 
will ignore this complication throughout this paper. 

The Abelian case is not entirely trivial, since the action S2 contains the square root of an operator — the hallmark of 
the FSE. (Throughout this paper, we take this operator, called f2, in the simple form £1 = \/M 2 — V 2 .) Nonetheless, 
S2 has center-vortex solutions. Although these do not completely coincide with conventional d — 3 center vortices, 
they show the necessary features: Long-range pure-gauge parts that confine, and field strengths that vanish at large 
distances as exp(— Mp) where p is the distance from the closed string on which the vortex lives. 

In anticipation of what we must do in the non- Abelian case, we study briefly the infrared expansion of £2 in powers 
of k 2 /M 2 , and show that the first two terms yield a familiar action. The leading term is a gauge-invariant mass term; 
the next-leading term is the usual Abelian gauge action. However, if the expansion is truncated after two terms, the 
gauge mass described by them is erroneous. The reason is elementary: The infrared expansion, at least in the Abelian 
case, is nothing but the first two terms of the expansion 

VM 2 ~V 2 i(M 2 - l -V 2 + . . .); (1) 

the two terms saved correspond to a mass y/2M instead of M. We propose that it may be phenomenologically useful, 
although not highly accurate, to make the replacement 

VM 2 -V 2 — -^(M 2 -V 2 ) (2) 

where the renormalization constant Z ~ 1 can be estimated in various ways. This heuristic replacement has the 
correct gluon mass. We discuss the motivation for this renormalization, coming from omitted terms in the infrared 
expansion. 

The second step, the subject of Sec. lIIIl begins with the problems of enforcing non- Abelian gauge invariance. Using 
earlier work 0], we show that S2 of step one can be gauge-completed to exact non- Abelian gauge invariance with an 
infinite series of n-point functions and powers of gauge potentials, in such a way that all n-point functions depend 
only on the operator O, no matter what the specific form of fl is. Gauge completion uses the pinch technique (la . [23| 
and the gauge technique [24j, as reviewed in the Appendix. The gauge technique is an approximation that becomes 
exact only at zero momentum, but is useful generally for momenta not large compared to M. In this gauge-completed 
S we continue to use the simple form in S2 of the two-point function introduced in the Abelian case. Just as for 
the ordinary Schrodinger equation, either direct substitution in the FSE or a dressed-loop expansion based on [25j 
ultimately yields a non-linear Schwinger-Dyson equation for fl (see the Appendix). We do not attempt to carry 
out this difficult program to find f2, but simply use the form already introduced in the Abelian case, showing mass 
generation. 

Even the approximate (although showing full non- Abelian gauge invariance) form of S coming from application of 
the gauge/pinch technique is extremely complex, involving not only square roots of operators but an infinity of terms. 
This S looks nothing like actions that we are used to dealing with. Ultimately, whatever form S takes must be dealt 
with on its own terms. However, just as in the Abelian case it can be helpful to look for an approximate but familiar 
form. We make the same sort of mass expansion, saving only the first two terms, and argue that for QCD the leading 
term, of 0(M), is equivalent to a gauged non- linear sigma (GNLS) model, which is commonly used as a description 
of gauge-invariant dynamical mass generation in Yang-Mills theory (see, for example, (l8l . [2q|). This sigma model 
contains the massless scalar poles, actually pure-gauge parts of center vortices, that are responsible for confinement. 
The second term, of 0(1/M), is (after gauge completion) the conventional Yang-Mills action. But as in the Abelian 
case, the mass is wrong by a factor y/2, so we suggest using the replacement of Eq. ©. 

In Sec. II VI we give the final conjecture for the non- Abelian exponent S and d — 3 action Id—s = 2S, and the main 
consequences following from it. The conjectured action is the sum of a GNLS model and a conventional Yang-Mills 
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action, with the correct free-field mass and a poorly-known renormalization constant Z . We suggest a method or two 
for estimating Z , probably with no more than 25% accuracy. 

The fourth step is to examine the consequences of this final two-term action. We have already noted that this 
action has center vortices as classical local minima (classical maxima of the FSE wave functional), and thus could 
provide a description of confinement in the FSE, which was one of our principal goals. Moreover, by appealing to 
known d = 3 gauge dynamics, we can estimate the d — 4 coupling strength in terms of the renormalization constant 
Z . In d = 3 the coupling g\ has dimensions of mass, and there is a unique [f or g iven SU(N) ] dy namically-determined 
ratio M/g%, which has been estimated by a number of authors [13, [H, [29|, Ho, HH, [H, y!, Hj, HH, UlI Knowing 
only this ratio we can estimate the d = 4 QCD coup ling a s (M 2 ), getting a value around OAZ. For Z ~ 1 this is 
reasonably close both to an early d = 4 estimate [18[ using the gauge technique and pinch technique but somewhat 
low compared to phenomenological estimates [2l[ of 0.7±0.3. Another application is to the d — 2 + 1 FSE, studied in, 
among other works, @,[I3|- Our present techniques suggest that the corresponding d = 2 FSE exponent S is again 
a sum of a gauge-invariant mass term and the usual Yang-Mills action. Greensite speculated that this S just had 
the conventional Yang-Mills term. However, as noted there and in [T3], this would lead to the wrong conclusion that 
in d = 2 + 1 all representations of SU(N) were confined, when in fact the adjoint and other representations with 
iV-ality = mod N are screened, not confined. But with the addition of the mass term, confinement can come about 
through center vortices, and this form of confinement correctly predicts screening for these representations. 

The paper ends with Sec. [Vj giving conclusions. An Appendix reviews some background material on the FSE, 
including applications of the pinch/gauge technique to the gauge FSE. 



II. DESCRIBING MASS GENERATION IN THE FSE: THE ABELIAN CASE 

Notation: Throughout this paper we will always use the canonical gauge potential Af(x) potential multiplied by 
the coupling g, with the notation: 

At(x) = gAm. (3) 

Here a is a group index for gauge group SU(N), and i — 1, 2, 3 index the spatial components. All vectors are three- 
dimensional, so we will now drop the vector notation and just use, e.g., k for a three-momentum. We also use the 
antihermitean matrix form 

M*) = (£)AaA°(s) ( 4 ) 
where the A a are the Gell-Mann matrices for SU(N), obeying 

Tr \ Xa \ Xb = \ 6ab ' ^ 

The A°~ have engineering mass dimension 1 in any dimension. The time component A® is missing from the FSE. In 
this paper we will not need to indicate gauge-fixing and ghost terms necessary to define the d = 3 functional integrals 
that yield physical expectation values. 

In the first step we begin with a simple quadratic (in the gauge potentials) form for S that is consistent with gluon 
mass generation. This quadratic form S2 is Abelian, showing U(1) N _1 local gauge invariance: 

1 

27 



S i = ^2 I A-l^A^x) (6) 



where the integral is over three-space, and Qij is a product of two factors: 

Q i j = p ij n. (7) 

The factor Pjj is a transverse projector: 

A, = - ^ (8) 
that is required for Abelian gauge invariance. The free-field value of f2, called Oq, describes free massless particles: 



fio = V 3 ^ = VP (9) 
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where k is a three-momentum. To describe dynamical mass generation we will use, in this paper, the simple form 

CI = V-V 2 + M 2 (10) 
in which the gluon mass M is the on-shell value of a running mass. Putting these equations together we have: 

S2 = WI A i^ M2 -^ A J- C 11 ) 

One can easily check that S2 is an exact solution to the FSE for an Abelian Hamiltonian with a gauge-invariant 
mass term: 

H = l{-lg 2 (^) 2 + ^[l(Ft 3 ) 2 + M 2 AtP t] A^}^ j[\w) 2 ] + V. (12) 

where = diA^ — djAf are the Abelian field strengths. Here the mass term is put in by hand; in the non-Abelian 
version, we imagine that this mass term summarizes the effects of non-Abelian condensates. 



A. Equations of motion and solitons for S2 

One goal in this Abelian example is to find center vortex- like solitons as extrema of 5*2. It may not be entirely 
obvious how to proceed, because this action has the square root of an operator, leads to s ubtleties concerning positivity, 
locality, and self-adjointness. For example, we will see that the operator V M 2 — V 2 effectively vanishes on center 
vortex solitons, although —V 2 is formally positive; this would falsely suggest that the action of such a soliton is zero. 
Consider the following alternative description of S%, found by expanding the square root in powers of — V 2 /M 2 and 
assuming that integration by parts with no boundary terms is allowed at all orders: 

S2 = ^J AtPijA i + 4? / £ Cn+iM- 1 - 2 "^ . . . d N F«{x)] 2 (13) 

where dk = d/dxk and the CV are the coefficients of x N in the power-scries expansion of VI + x - This re-definition 
of the square root gives the same generalized Euler-Lagrange equations as the naive equations following from the 
original form of Eqs. (I6|7I10|) . because these equations assume that integrating by parts gives no contributions (as 
would be appropriate for functions that fall off at least exponentially). 

In order to study these generalized Euler-Lagrange equations, it is very helpful to have 62 in a formally local 
form. We note that, term by term, all but the first term of this alternative form of S2 are both local and manifestly 
gauge-invariant, and need no change. As for the first term, we replace (in a familiar way) the non-local part by scalar 
fields: 

Mr 1 r 

S 2 = ^/[^-^ + ^/W + ... (14) 

Now keeping only a finite number of terms in the mass expansion of S2 yields a local action, although of course the 
infinite sum may introduce non-localities. 

Saving only the first two terms in the mass expansion of S2 based on Eq. (fT"3f should fail to satisfy the Abelian 
FSE based on the Hamiltonian of Eq. (|12jl . It is instructive to work out this failure and its consequences. The FSE 
reads: 



2 J "SAf 2 J SAfSAf 

where E is the vacuum energy. Since the second-derivative term on the left-hand side of this equation only contributes 
to E, we drop it and renormalize E to zero. The mass expansion of S2 suggests that the remaining quadratic term in 
the FSE is in error at C(l/M 2 ). A simple calculation confirms this; Eq. (fl~5|) becomes: 

^/<H> 2 + " + ^/2>^> 2 = 0. W 

At least qualitatively this error term in the FSE [last term on the left-hand side, of 0(1/M 2 ), the same relative 
order as the N = 1 term in Eq. (|13p] can be thought of as increasing the kinetic field-strength term (F^) 2 by a factor 
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involving a mean-square momentum of the type (k 2 )/M 2 ; such an increase helps restore the balance between kinetic 
and mass terms in the expanded Hamiltonian which was disrupted by the usual infrared expansion of Eq. (|14[) . Such a 
renormalization is not quantitatively trivial, since momenta relevant for solitons such as center vortices are of 0(M). 

It is useful to restate the local form of 52 in a compact way, by undoing the power-series expansion and integration 
by parts: 



S2 = wl [4 ° " dir]2 + wl A ' Pl3 [y/M * v2 " M]A >' 



(17) 



The scalar fields </> a are to be integrated over, which may be thought of as projection of a simple mass term (Af ) 2 
onto its gauge-invariant part by integrating over all gauge transformations. Because the <p a appear quadratically, such 
an integration is the same as solving the classical field equations. The field equations for the 4> a are identical with a 
constraint following from the field equations for the Af. 

Varying S2, one finds the gauge potential equations of motion: 



M (Af - d^ a ) + [VM 2 - V 2 - M I',, A; = 0. (18) 
The divergence yields the 4> a equations: 

1 

V 2 



V V = diAf -> <p a = =^d l Af + ip a with V V = 0. (19) 



Re-write Equation (TT9j) as: 



y/M 2 - V 2 P l3 Af = Md^ - =0^?) = Mdup*. (20) 



Multiplication by \J Af 2 — V 2 leads to: 



(M 2 -V 2 )PijA] = M\/M 2 - V 2 dnp a -> (21) 
V 2 Af-d i d j A^ = M 2 (Af - d~djA1) - M\[M 2 - V 2 d i <p a — 
W 2 Af - didjA? - M 2 (Af - d l( j> a ) = M[M - a/M 2 - V 2 ]d i cp a . 



Term by term, every term on the right-hand side of the third equation in Eq. (|2ip vanishes, if we use ~S/ 2 ip a = 0. 
Since in 1Z 3 there are no fields ip a solving V 2 y> a = that are regular everywhere and vanish at infinity, one may 
be tempted to make the stronger statement that ip a must vanish. But the description of center vortices requires 
a non-zero ip a , singular on a closed Dirac hypersurface of co-dimension 2 (a closed string in d — 3), so it is more 
accurate to say that term by term the expansion of the right-hand side of the third equation in Eq. (|21j) vanishes 
almost everywhere. However, we will soon see that this is not true for the unexpanded form. If we nonetheless drop 
this term with the square-root operator, the final equations of motion are the usual equations 38] for center vortices, 
the same as would be gotten from the d = 3 Euclidean action 

^ J {M 2 (A? - d t r? + \{F°-) 2 }. (22) 



This action is just the potential V occurring in the Abelian Hamiltonian of Eq. (|12| . written in local form; it is the 
Abelian form of the d = 3 infrared-effective action used [H, [38| to describe mass generation, and it has center vortices 
as classical solitons. 

If the term M[M — V M 2 — V 2 ]diip a is left unexpanded, things are slightly different, although there are still center 
vortices characterized by long-range pure-gauge parts and field strengths vanishing exponentially as exp(— Mp), where 
p is the distance from the Dirac string. A center vortex is always fully determined by tp a . We present our results in 
the gauge diAf = 0, in which case 4> a — (p a . The well-known expression [38( for the center-vortex di4> a is: 

d^ a (x) = 2nQ a e llk d 1 j> dz k J ^^A e *(-*) (23) 

where the closed contour T is the Dirac string, and Q a is one of the N — 1 generators of the Cartan subalgebra, 
normalized so that exp(27riQ) is in the center of SU(N). Now the third equation in Eq. (|2"T|) easily gives: 

r f d 3 k M 

Af(x) = 27rQ a e i7 - fe a- d> dx k / -— r e iHx-z) / 24 ) 

iy ' 3 7r J (2^) 3 k 2 Vk 2 + M 2 V ' 
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In the usual d = 3 vortex, an extremum of the action in Eq. (|22)l . the factor M(k 2 + M 2 ) x l 2 would be replaced by 
M 2 (k 2 + M 2 )- 1 . 

This unusual square root does not change the fact that the field strengths show exponential decrease; in fact: 



1 f M 2 

Bt = = 2vQ a j^ dZi _ Kx(M\x - z\). (25) 



}.f2 

There is, of course, still the long-range pure-gauge part associated with cf> a , which we can isolate by the decomposition: 



M _ 1 1 M 
k 2 s/k 2 + M 2 ~k 2 + k 2 Wk 2 + M 2 ~ 1 *~ (26) 

The second term on the right-hand side is short-ranged. The short distance behavior is more singular than that of the 
conventional vortex, but leads only to a logarithmic singularity in the value of S2, the same as for the conventional 
vortex. In both cases the singularity is multiplied by a power of M , which removes the singularity because the running 
mass vanishes at short distances. So the vortex extrema of S2 differ in detail from the usual center vortex, but have 
the hallmark features of a long-range pure-gauge part and field strengths vanishing like exp(— Mp). 



B. Mass expansion of S2 

Another goal of this section is to replace S2, which is either given in Eq. (fT3|) as an infinite sum involving derivatives 
of arbitrarily high order or in Eq. (1171) in terms of square roots of operators, by a tractable and recognizable action. 
The first two terms of the expansion, written explicitly in Eq. (|14l fit these criteria, but suffer from a serious defect. 
The coefficient of the second term, the usual gauge action, is wrong by a factor of 2; as written, it describes gauge 
bosons of mass \/2M. This wrong coefficient arises from the expansion \/l + x = 1 + (x/2) + . . .. We can see the 
same thing happening with a mass expansion of the Fourier kernel of Eq. (|2l)|) . Expand the square root in the curly 
brackets of this equation in powers of k 2 /M 2 to get: 

1 1 (27) 



fcVfc 2 + M 2 k 2 k 2 + 2M 2 

This is exactly the kernel of the usual d = 3 vortex, but with the wrong mass \J~2M . This is not the only way of 
expanding; for example, re-writing the Fourier kernel in a different form and expanding the square root occurring in 
it gives: 



M MVk 2 + M 2 M 2 1 

(28) 



k 2 Vk 2 + M 2 k 2 (k 2 + M 2 ) k 2 {k 2 + M 2 ) 2(k 2 + M 2 ) 

The first term on the right-hand side is the standard center vortex with the correct mass M, and all other terms 
have this mass as well. However, these other terms give the wrong coefficient for the exponential falloff of the field 
strengths at large distance. 

There are no such results for square-root operators in the non-Abelian case, which is as expected much more 
complicated. So we will, in the spirit of the Abelian expansion given in Eq. (TH| . look for a way to approximate the 
complicated non-Abelian result by a two-term form, the first of which is a (gauge-invariant) mass term and the second 
is the usual Yang-Mills action. In the Abelian case, such a two-term action as an approximation to the infinite sum 
of Eq. (fT"3"|) suggests that the derivatives in this sum, beyond those in F°- , be approximated by averages so that this 
equation is effectively 

MP '\ f k 2N M f 7 f 

52 = 2i 2 I AtPijA 1 + W J E C ^^M^ )[ ^ (X)]2 S V J + 4M? J m 

where k 2N stands for the multiple derivatives. If this is justified, the infinitely many terms of Eq. (|13p are indeed 
replaceable by a mass term plus a renormalized conventional gauge action. But because the gluonic mass described 
by this S2 must be M, the same as in the original S2, there will have to be an equal renormalization of the mass term 
in Eq. (|29[) above. In later sections we will explore an approximation to the square root that is motivated by these 
remarks, involving the replacement 

Z 



^/M 2 - V 2 -» jj(M 2 -V 2 ) (30) 
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for some renormalization constant Z, supposed to be near unity. We have no reliable techniques for calculating Z, 
so we will resort to a simplistic approach of making a least-squares fit of the operator \JM 2 — V 2 by the operator 
(Z/M)(M 2 - V 2 ), which leads to Z ~ 1.1 - 1.2. 

Before engaging in this mass expansion we must understand the gauge structure of the non-Abelian exponent S. 



III. THE NON-ABELIAN CASE: GAUGE COMPLETION AND MASS EXPANSION 



We can be nowhere near as complete in the non-Abelian case as we were above, and ultimately will be forced to 
resort to a large-mass expansion. 

In the non-Abelian case, the quadratic term S2 with which we began is supplemented with an infinity of terms, 
involving spatial integrals over n > 3 spatial gauge potentials multiplied by an n-point function f2„ depending on the 
spatial and discrete coordinates of the gauge potentials (see the Appendix): 

g2s = kJ I m ^ + h.j j j + ■ ■ ■ (31) 

The n-point function of this expansion is related to the n + 1-point function through ghost-free Ward identities, as 
arise in the pinch technique [H, |2j|. These Ward identities can be "solved" using the gauge technique, a well-known 
technique whose main points of interest we describe in the Appendix, and the result is that it is possible in principle 
to find an approximate but exactly gauge-invariant expression for the entire series of n-point functions describing the 
wave functional exponent S. Each n-point function depends only on the two-point function, but in a complicated way 
that is not understood. Ultimately the two-point function is determined by a non-linear Schwinger-Dyson equation 
that can (again in principle) be derived either by direct substitution in the FSE or by a dressed-loop expansion 
[25l [39l . |40 |. Using a dressed-loop expansion for S is equivalent to direct solution of the FSE (of course, either the 
dressed-loop expansion or the FSE must be truncated at at certain number of loops, but this truncation has nothing 
to do with a truncation in the coupling g 2 ; all-order non-perturbative effects arise even at one-dressed- loop order in 
QCD). 

A systematic study of the FSE would go on to determine the mass M from the infinity of equations for the n-point 
functions in S of Eq. (|3"Tj) . but that is not our purpose here. Instead, we show how to construct what we will call 
a gauge completion of the 2-point action for an arbitrary f2, using earlier work Q, to add higher-point functions, 
consistent with solving the FSE, that depend on f2 in specific ways that insure full non-Abelian gauge invariance. 
Ultimately, the FSE becomes a non-linear equation for f2, just as for the ordinary Schrodinger equation. 

Gauge invariance requires that the lowest-order (quadratic) term has the Abelian form already given in Eq. pip . 
The Ward identities for the three-point function and their solution are detailed in the Appendix. Both the Ward 
identities and the FSE for the determination of this three-point function involve only the two-point function fly, 
and it is plausible that there exists a three-point function satisfying these equations that is a functional solely of the 
two-point function fly. The gauge technique provides such a three-point function, as given in Eqs. (|A20IA2ip . The 
gauge technique by itself does not furnish a unique solution, which must be found by recourse either to the FSE itself 
or to the dressed-loop expansion. However, in the infrared limit of momenta small compared to the mass M the 
solution is unique. 



A. Mass expansion: The leading term 

In general, the gauge/pinch technique leads to quite complicated expressions, and we will explore only a simplified 
version of it. The main simplification is to look at the leading terms in an expansion in inverse powers of M. In the 
leading term, of 0(M), all two-point functions f2 are replaced just by M itself, which gets rid of many momentum- 
dependent terms. In this way the leading term of the three-point function is: 



^HhJ^ks) = f abc Yi kuk2 %i k2)k + '■•"■> + wv.. (.->>->! 



One can proceed in principle this way, by looking at the pinch/gauge technique solution for the four-point function 



(see |4l|) and taking the large-mass limit, then the five-point function, etc. We will not detail such an investigation 
here, but will point out some features that strongly suggest the all-order solution. The structure of the Ward identities 
shows that the leading term of every n-point function is 0{M), with all other dimensions taken up by momenta, and 
that the gauge-technique solution involves longitudinally-coupled massless poles whose number grows with n. Observe 
further that the GNLS term of S is the exact solution of an FSE Hamiltonian consisting of just this term itself, as 
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given in Eq. (|33[) below, multiplied by M. Of course, there is no such term in the underlying QCD Hamiltonian, but 
there would be one in the infrared-effective Hamiltonian of QCD, derived by d = 4 techniques [l8l l26l|. 

We suggest that the action of the gauged non-linear GNLS model, expressed in non-local form [as in the originally- 
stated form of 5*2, in Eq. (TTTj) ] is the all-order perturbative solution to the leading mass terms of the gauge/pinch 
technique approach. To find this non-local form we investigate the classical solutions of the local GNLS action. 

Because the notation is more compact, we temporarily switch to the antihermitcan matrix notation of Eq. (j4|). The 
local GNLS model, normalized appropriately, has the action [26j | : 

— M r 

Ignls = -jjTj ^xTriU^DtU} 2 (33) 

where U is a unitary matrix transforming as U — > VU under the gauge transformation 

Ai -> VAiV' 1 + VdiV~ x . (34) 
The classical equations for U express this quantity in terms of the Ai [IB], with the result 

U = e»; u=^d.A+±{[A i ,d i ±d.A] + ±ld.A,±d-A]+:-} (35) 

showing the appearance of massless scalars. More generally, since U~ 1 DiU is a gauge transformation of Ai, functional 
integration over the U is equivalent to projecting out the gauge- invariant part of the mass term [1, [ll|. Note that the 
term linear in Ai of the GNLS model field U~ 1 DiU is the transverse part of A». This linear term is Abelian, and all 
higher-order terms of u in Eq. (|35j) are non- Abelian. 

[Greensite and Olejnik [TtJ have conjectured that in certain instances operators such as V~ 2 should be replaced 
by D~ 2 , where = di + Ai is the covariant derivative. Their lattice calculations show that D~ 2 is a finite-range 
operator, with no massless poles; this is reasonable, because it contains gauge-potential condensate terms, but it is 
not obvious where the long-range pure-gauge excitations responsible for confinement, such as we have in Eq. (I35[) . 
are. We will not follow this line of reasoning here.] 

It is now straightforward, if tedious, to verify that the two- and three-point terms of the non-local GNLS action 
give rise precisely to (the leading mass terms of) the two-point function 5*2 and the pinch/gauge technique three-point 
function of Eq. {32]). Moreover, the GNLS action integrated over U automatically satisfies the Ward identities to all 
orders, just because it is the solution of the FSE for a gauge-invariant Hamiltonian. 



B. The second-leading term 

We already know that the next-leading term, of 0(1/M), in the expansion of the two-point function S% is the 
conventional Abelian action involving Ffj. It is obvious without any calculation that the Abelian action will, at a 
minimum, be gauge-completed to the full Yang-Mills action with its three- and four-point vertices. These come from 
the three- and four-point functions in the expansion of S as given in Eq. (|3"Tj) . The desired terms of the Yang-Mills 
action are straightforwardly found either by direct solution of the FSE or from the dressed-loop expansion, which 
always contain all the terms of the action of the underlying theory divided by some sum of two-point functions O. 
For example, we show in the Appendix that the three-point function has the term 

Qflgfa, k 2 , k 3 ) = [0(1) + 0(2) + 0(3)]- 1 / abc fe(A;i - k 2 ) k + c.p] + ... (36) 

where the term in square brackets is the free Yang-Mills three-point vertex and each 0(i) is replaced by M to find the 
leading term in the mass expansion. There is a plethora of other terms, which either cancel among themselves or give 
total divergences. Of course, higher-order gauge-invariant terms may arise from higher-order coefficient functions in 
the gauge-potential expansion of S, Eq. (|A14|) , but we will not consider them, since they are necessarily accompanied 
by higher powers of 1/M. 

In the Abelian case the 0(1/M) term is of the correct functional form, but with a coefficient twice as small as it 
should be, and the same problem arises for the non-Abelian case. This results in a gauge mass of \/2M instead of M, 
as pointed out in Sec. [Til I n the next section we consider a modification of the straightforward mass expansion of the 
type of Eq. (|3T))) that forces the correct mass. 
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IV. THE FINAL CONJECTURE AND ITS CONSEQUENCES 
A. Heuristic mass expansion 

What we have so far in the gauge-completed mass expansion to second order is the sum of a GNLS and a Yang-Mills 
term, but with the wrong mass. What we need is an approximation to this two-term action that has the correct mass, 
in part because solitons are described in this momentum range and decay at a rate ~ exp(— Mp). In any event, it 
is clear that the first two terms in any sensible infrared expansion consist first of a gauge-invariant mass term and 
second of a standard Yang-Mills action. 

Rather than stick to a strict expansion in powers of V 2 /M 2 , we conjecture that, as in the Abelian case, we can 
replace ft = V- V 2 + M 2 by a leading term (Z/M)(-V 2 + M 2 ), where Z is a coefficient of 0(1). 

The mathematical motivation for least-square fits of operators is well-known. Consider a normal operator P, 
expressed in terms of its eigenvalues and eigenfunctions: 

P = £>>*nR (37) 

Any function of P, call it f(P), is expressed by replacing X n by f(X n ). With the operator norm TrP^P, we define a 
relative RMS distance between two operators f(P) and g(P) by: 

/ Tr[f{P) g{P)] [f(P) - 5 (P)]t | 1/2 _ J / dAp(A)|/(A) - g(X) | 2 \ 1/2 ^ 



I Trf(P)f(Py J \ JdXp(X)\f(X)\ 2 

where 

P (X)=J2S(X-X n ) (39) 

is the density of eigenvalues. One could also modify this density by multiplying it by a non-negative function q(X) to 
emphasize a certain range of eigenvalues, so that the weight in the integral is p(X)q(X). 

The eigenvalues of P = —V 2 are the squared momenta fc 2 , positive for real fc. We really want our approximation 
of ft to be fairly good for imaginary fc, so the above discussion is not very useful. Moreover, the operators involved 
are not in trace class, so divergences arise. Instead, we take a rather simpleminded point of view, asking what is 
the best fit, in the least-squares sense, of the function Z(l — x 2 ) to the function yl — x 2 in the interval < x < 1. 
Here x 2 represents V 2 /M 2 , and positive values for this operator suggest that we are applying it to a special class of 
functions representable by Laplace transformation, with the Laplace-transformation weight peaked around M. This 
is indeed the property of the functions that enter into FSE center vortices, as exemplified in the Abelian center vortex 
of Eq. ([25]). 

For a uniform weight over < x < 1 we find the normalized least-squares integral Ii x (Z): 



hs(Z) 



^dx[Z(l-x 2 )-^/T^x 
Jo dx Q - x2 ) 



-[ 1/2 

(40) 



Minimizing on Z gives Z = j^l ~ 1.10, and the minimum value of Ii s is about 0.22. If we replace x 2 by x in the 
integrand of Eq. (|4"0")) . which corresponds to a different weight, we get Z = 1.2. Both values are near unity, as expected, 
and the value I/ s ~ 0.22 suggests the relative accuracy of this least-squares fit. 

B. The final conjecture: Relating d- and d — 1-dimensional dynamics 

The final form of the conjecture, expressed in terms of the d = 4 variables g 2 , M is then: 

-2S= -J d=3 - ^ J d'xTriU-'D.U} 2 } + f d?xTrG% + 0(M" 3 ). (41) 
We now compare this to the canonical d — 3 form of the conjectured action, action, which is: 

/d=3 = - / d^xi^TrGl + ^Tr[U- l DiuA . (42) 
J I z 9s 93 J 
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Here Gij is the non-Abelian field strength, Di — di + Ai is the covariant derivative, and the unitary matrix U is the 
GNLS field, as before; the gluon mass is 7713 and the d — 3 coupling, with dimensions of mass, is gf. Equating Id=3 
with 25 leads to: 



ZM _m\ 2 _ 2Z 93 
These equations yield m 3 — M, as expected, plus 



(43) 



9 2 = (44) 



TO 3 



(Note that the d = 3 quantities scale properly at large TV if their d — 4 counterparts do.) Presumably the d = 4 
coupling g 2 that occurs in these formulas is actually the running coupling g 2 (M 2 ) evaluated at the gluon mass scale. 

We can now make an estimate of a pure d = 4 quantity in terms of a pure d = 3 quantity, from Eq. (|44[) and earlier 
d = 3 results. In d = 3 one quantity of particular interest is the dimensionless ratio m^/g 2 . This ratio has been 
estimated in a number of continuum and lattice studies [U H3> [2t| [3(| HH, HU, [H, [13, [H[ , and we can see whether 
this d = 3 dynamical quantity can correctly predict the running coupling g 2 (M 2 ) at the gluon mass scale. Or we can 
reverse the problem and use estimates of the running coupling to predict m^/g 2 . There is no particular reason to 
think that the dynamics of the action defined by the exact vacuum wave functional, before truncation to two terms 
of a mass expansion, should be precisely that of d = 3 QCD. Nonetheless, if our conjecture is to be believed there 
should not be gross discrepancies. 

In SU(2) gauge theory various authors [13, [H, [2t| [3(| HH [H, [H, [H], [H[ give a value m^/g 2 ~ 0.32, and one 
SU(S) lattice study [36[ gives a value of 0.48. The quantity ma/gf should be linear in N of SU(N) for large N, 
and the factor 3/2 nicely converts the 517(2) values to the SU(3) value, so we use 0.48 as the SU{3) value. We 
then find a value for the strong coupling (with no quarks) a s (M ) ~ 0.33Z that is in fairly good agreement with 
the one-dressed-loop approximation found in the original paper on dynamical gluon mass generation [181 ] . This paper 
gives a one-dressed-loop equation for the running charge with dynamical gluon mass generation. At the momentum 
scale of the gluon mass M: 

a ^ = i = m^^MW)] * °- 4 ' (45) 

where the numerical value is based on the estimates M — 0.6 GeV, A = 0.3 GeV, and the absence of quarks (Nf = 0). 
Of course, these numbers for M and A are themselves uncertain, if only because Eq. (|45|) is a one-dressed-loop 
equation. 

According to this one-dressed-loop equation, accounting for three light flavors multiplies the no-quark value by 
11/9 ~ 1.2. If we assume that this correction applies to the FSE result of this paper, which as it stands does not 
account for quarks, our estimate of a s (M 2 ) increases to about 0AZ. 

Several papers have extracted values of a s (0) ~ 0.7 ± 0.3 from various scattering data sensitive to low-momentum 



effects [2l| that could diverge if there were no gluon mass. The three-quark value that we give of 0AZ is a little 
smaller, but in quite reasonable agreement considering the approximation that is inherent in a two-term truncation 
of the FSE exponent 5 and our lack of knowledge of Z.. 

It has been argued [37] that m^/g^ for SU(N) is very closely approximated by the simple analytic function 

^ = — ■ (46) 

the present author [ll| has argued for a ratio that should be fairly close to 15iV/(327r), which differs from the above 
by only a few percent. One then has a simple analytic formula for a s (0). Using the value from Eq. (|46|) in Eq. (|44|) 
yields the amusing, if not very accurate, formula 

aJM 2 ) = — -—. (47) 

We can play the same game in one less dimension for the d = 2 + 1 FSE, beginning with an exponent 5 for the 
wave functional that is the trivial dimensional reduction of what we began with in d = 3 + 1. The result is a d = 2 
action with, as in d = 3 + 1, a mass term and a kinetic term. This is not the standard d = 2 QCD action, which is 
a free field theory. We compare this to a conjecture made long ago by Greensite 0, arguing that 5 for the FSE was 
just the usual Yang-Mills action in one less dimension. Unfortunately, as [13] notes, if Greensite's 1979 conjecture is 
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applied in d = 2 + 1, the effective action is the familiar d = 2 free-field QCD, which would lead to confinement of 
all representations of SU(N), not just those with Y-ality nonzero. This is not the right behavior for d = 2 + 1. But 
in our case once again the action is the Yang-Mills term plus a GNLS model mass term; this action has [3§| center 
vortices; they are point-like objects in d = 2. A condensate of these solitons leads to confinement, but only of group 
representations that have Y-ality ^ mod N; other representations (such as the adjoint) are blind to the long-range 
parts of center vortices. This is the correct behavior for d = 2 + 1 gauge theories. However, if the mass term were 
not present in ld=i this action, which is supposed to carry all the information about d = 2 + 1 gauge theories, would 
reduce to the standard Yang-Mills action in d = 2. The conventional treatment of d = 2 gauge theories, which (in 
the absence of dynamical matter fields) are free field theories, finds confinement through the long-range free gluon 
propagator, and all representations are confined. But with the mass term the gluon propagator is short-ranged and 
confinement comes from the pure gauge long-range parts of center vortices. 

It is far from trivial to calculate the properties of the center- vortex condensate in d = 2, and so we cannot relate 
the d = 3 coupling to the string tension that would be found from the d = 2 effective action. 

V. SUMMARY AND CONCLUSIONS 

We have conjectured that to a reasonable approximation the dominant quasi-infrared part of the vacuum wave 
functionals for the d = 3 + 1 and d = 2 + 1 FSE are actions in one less dimension consisting of a Yang-Mills term and 
a GNLS model term, showing gauge-invariant dynamical mass generation. Two main conclusions follow: 

1. Given the usual entropy-dominance argument, these wave functionals show confinement through center vortices, 
such that only group representations with Y-ality ^ mod N are confined. 

2. In d = 3 + 1 we can appeal to earlier works estimating the ratio 7713/ '<?f in the d — 3 action of the FSE to make the 
estimate a s (M 2 ) ~ OAZ, where Z is a renormalization constant that we have very crudely estimated to be in the 
neighborhood of 1.1-1.2. This can be compared to an earlier estimate, based on the original work on dynamical 
gluon mass generation, of a s (M 2 ) ~ 0.4. Both these estimates have three light flavors of quarks. This is to be 
compared to compared to phenomenological estimates [2l|, also with three light quarks, of a s (0) ~ 0.7 ± 0.3. 

It would be interesting to verify this structure of the FSE vacuum wave functionals through lattice simulations. 
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APPENDIX A: A BRIEF REVIEW OF THE FSE AND SOLUTION METHODS 

The purpose of this review of known material Q is to indicate the plausibility of constructing an infrared-accurate 
and gauge-invariant form of the wave functional -0, based on a single operator ft, obeying a non-linear Schwinger- 
Dyson equation. In ordinary quantum mechanics this is exactly what happens except that the "Schwinger-Dyson 
equation" is simply algebraic. The FSE for scalar field theories is really nothing but ordinary quantum mechanics for 
infinitely-many coupled oscillators, so we review it and its connection to quantum mechanics, then go on to gauge 
theories. 

1. The Schrodinger equation 

The general principles of solving the FSE in terms of an operator ft are most easily understood from the ordinary 
Schrodinger equation. Consider the quadratic/quartic Hamiltonian 

The ground-state solution is ip = exp(— S), with 

S= ^Qx 2 + ^Qix* + . . . (A2) 
Following Q, we substitute ip in the Schrodinger equation saving only terms through ftg and find: 

-5fi| A Oe n a 2 + i n E = 1 (A3) 

6 3ft ' 4 4ft 8ft' 4 4 ' 2 V ; 

It is easy to solve the equation for ft4 to derive a quartic equation for ft. One can go on to any order this way, 
expressing (in principle, at least) every n-point coefficient up to a given highest value of n in terms of ft, and ending 
up with a non-linear dressed-loop equation for ft. Consider now the case uj — 0, for which the perturbative expansion 
coefficient X/oj 3 diverges. Then through the six-point term we find the expression E = (3A/272) 1 / 3 , which has the 
value 0.2226A 1 / 3 . This is within a few percent of the numerical answer of 0.2311A 1 / 3 . 

2. Field theories other than gauge theories 

For simplicity of exposition, we begin with a scalar field theory. Take the FSE Hamiltonian to be 



H = d 



j3 . 



(A4) 



where the potential V contains cubic and higher terms. Ref. [8j showed that the vacuum wave functional could be 
expressed as a d = 4 partition function: 

e~ s = const, x /(dfc) exp[-/ ($) - I ($o) - J V{<P + fo)]. (A5) 

In this partition function, space-time integrals are of the form of an integral over a Euclidean time t and all of three 
space: 

r°° r 

dr / d 3 x. (A6) 



o J 

The argument of S is the field 4>(x), and the field <t>{x) depends oni= (r, x) as: 

Mx) = e- n ° T <f>(x) (A7) 

with 



fto = JM 2 - V 2 . (A8) 
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The free action Iq is: 

Io($) = l |p T $)2 + (v$) 2 ] (A9) 
and the inverse of the free-action operator is the free propagator 

A = (z|-^-[e- n °l T - T 'l - e- n ^+^]\x'). (A10) 
The first term in the propagator is the usual Euclidean propagator: 



\x\——e 



1 , 1 f p ik-(x-x') 

~- n ° lT - Tl \x') = 7^ d^TT—T- (AH) 



2fi 1 ' (2tt) 4 J k 2 

For purposes of calculating the energy eigenvalue, this is the only term that needs to be saved in Ao, but the second 
term of the propagator in Eq. (|A10[) is necessary for calculating the wave functional. 

Either by working out the partition function of Eq. (|A5D or by direct substitution in the Schrodinger equation one 
sees that the vacuum functional ip has the general form (using a streamlined but transparent notation): 



4> = e- s :, S=\ [ U^ + Y.JTl [■ ■ /"<Wi...tor. (A12) 



2 J J ^— ' N\ 

j j N 

For purely three-dimensional equations, such as this, the unadorned integral sign simply indicates J d 3 x, where x 
is the argument of a corresponding (f> (and a sum over discrete indices, if any), with Q and the Qn,N > 3, as 
translationally-invariant form factors in the arguments of the </>. The partition function form in Eq. (|A5I c & n be 



addressed with the well-known resummation techniques 25j of the dressed-loop expansion. The effect of these rules 



is to remove a large fraction of one-particle-reducible graphs, as required for the dressed-loop expansion. In part, 
this amounts to a general replacement (but not quite everywhere) of the free operator f2o by a dressed operator f2 
that satisfies a non-linear Schwinger-Dyson equation. This operator is precisely the same as the f2 that occurs in the 
quadratic term of the wave functional in Eq. (|A12[) . 

For further details of this formalism for scalar field theories, see [4(| which uses it for calculating some terms in the 
Wigner distribution function. 



3. Gauge theories 



For gauge theories the same general structure holds; the principal problem remaining is to enforce gauge invariance. 
The canonical momentum and Hamiltonian are represented by 

n < - / + 27< B - > 2 i < ai3 » 

where Bf is the chromomagnetic field strength. The generator of infinitesimal gauge transformations is Dj b (— iS/SA b j), 
and this must annihilate -0. The exponent S in ip has the form given in Eq. (|3ip . repeated here for convenience: 

g2s = kj I • l < a < l < + h.j j I A< z Ah i A * n tk + ■■■ ( Al4 ) 

Invariance of S under infinitesimal gauge transformations is trivial for the two-point function fi^ ; this quantity must 
be conserved, so that in Fourier space 

n ij {k)=n(k)P ij (k); I',, .\, '-^f. (A15) 

For the free theory flo(k) = k. 

Gauge invariance is more complicated for higher-point functions. Annihilating tp with the generator of gauge 
transformations yields a set of ghost-free Ward identities (these Ward identities also apply to the pinch technique 
[~L8l [23| | construction of gauge- invariant Green's functions). For example, the Ward identity for the three-point 
function is: 

k u a&(ki, k 2 , k 3 ) = r bc [n Jk (2) - n jk (3)] (Aie) 
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where Q,j k (2) = £ljk(k 2 ), etc. 

Now turn to the FSE itself. The equation determining the three-point function has the general form 

n a mt% + n 3 -,( 2 ) n ** + ^(3)^ b = f abc r ijk . (An) 

The right-hand side Ty^ comes from the cubic term in H , plus another term from the five-point function. The Ward 
identity for r^j. is determined by the above equation plus the Ward identities for the two- and three-point functions 
as already given, and multiplying both sides of Eq. (|A17|) by kn yields: 

k u T ijk =n%(3)-n%(2). (A18) 

For free particles, with Q = f2o, this is satisfied by the usual free three-point vertex 

T° ijk = i(h - k 2 )kSij + c.p. (A19) 

The reader can verify that the FSE equation (1A17|) has a solution of the form: 

n#E(fci, k 2 , fc 3 ) = + n(2) + ft(3)r7 afcc |r yfc + {fi(i)^[n Jfc (2) - n ifc (3)] + C . P .}| (A20) 

which respects the Ward identity of Eq. (|A16|) . by virtue of the massless pole terms of Eq. (|A20I) . It should now be 
clear that these longitudinally-coupled massless excitations will occur, as a result of enforcing gauge invariance, for 
every n-point function. We will shortly identify these with couplings of the GNLS field introduced in our conjecture 
for the infrared-effective action. 

So far the vertex function Tjj k is undetermined. As [1| argues, one can carry out a program of expressing all higher- 
point functions in terms of the two-point function, and then the FSE (or the equivalent dressed-loop expansion) 
becomes a non-linear, non-perturbative equation for this two-point function CI. The idea, known also as the gauge 
technique, is to find an infrared-effective approximation to Tij k that exactly satisfies the Ward identity (|A18[) for any 
f2. One can, at least in principle, find such infrared-effective approximations for four- and higher-point functions as 
functionals of il. In fact, a very general form for the "solution" to the Ward identity for the three- and four-point 
functions is known [24], |4l[ for arbitrary dependence of ft on momentum. The word "solution" is enclosed in quotes 
because it is not unique; any completely-conserved term can be added to the "solution" for Tij k , for example. But 
the point is that purely-conserved terms are of higher order in momenta than the terms saved in the gauge technique. 

The general solution of |H) is: 

T m = Stjih - k 2 ) k - - fe) / n /fc (fc 3 ) - [P i i(k 1 )n lj (k 2 ) - Pjiik^Uuih)}^ + c.p. (A21) 

where the first term on the right is the free vertex r? fc and Ily(fc) = Pij(k)I\)k) is the transverse pinch-technique 
[HI, [23[ self-energy, related to fiy by: 

n^ = Py[ng + n{n}] (A22) 

where fi 2 , = k 2 is the free gluon contribution. 

In the simple case studied by us, II = M 2 , and the resulting expression for Tij k is: 

^ijk = dij{ki - k 2 ) k + — i-2-2 Yep. (A23) 

As we saw above, in ordinary quantum mechanics in one spatial dimension x the exponent S(x) of ip can be 
determined systematically from a set of non-linear algebraic equations, such that each term of 0(x 3 ) or higher can be 
expressed in terms of the quadratic coefficient tt. Finally, £1 is determined by a single non-linear equation, equivalent 
to a dressed-loop expansion. Combining the pinch technique and the gauge technique gives a completely analogous 
program for gauge theories, based on "solving" the Ward identities insuring gauge invariance. While this program can 
only be carried out approximately, it is gauge-invariant by construction. Ultimately it yields a dressed-loop equation 
for a single transverse operator f2jj(fe) = Pij(k)Vl(k). The pinch-technique self-energy II is itself a complicated 
function of f2, found by using dressed propagators of the general form given in Eq. (|A10[) . with f^o replaced by O 
and with appropriate vector kinematics. In effect, II is the on-shell self-energy and any ±iJ7(fc) occurring in If is a 
fourth component of a Euclidean four- vector (k 4 ,k) that is on-shell, by which we mean that k\ + k 2 + M 2 — 0, or 
fe 4 = ±Wk 2 + M 2 . 

All that we need from this development in the main text is Eq. (|A23[) . which will be used in the large-M expansion 
of the three-point function fi?^. 
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